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The statistical mechanics of a two dimensional Coulomb gas confined to one dimension is 
studied, wherein hard core particles move on a ring. Exact self-duality is shown for a version of the 
sine-Gordon model arising in this context, thereby locating the transition temperature exactly. We 
present asymptotically exact results for the correlations in the model and characterize the low and 
high temperature phases. Numerical simulations provide support to these renormalization group 
calculations. Connections with other interesting problems, the quantum Brownian motion of a 
QQ , particle in a periodic potential and impurity problems, are pointed out. 
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The physics of the inverse squared exchange Heisenberg model, the so called Haldane Shastry model [|l|j2| has been 
a field of considerable activity recently. Here one considers a spin i Heisenberg antiferromagnet in 1-dimension (1-d) 
on a ring of L sites with an exchange Hamiltonian 
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with Xj as integers denoting lattice points and = tt/L. The ground state of the model of Eq.(|l|) in a sector with N 
, spin reversals (relative to the ferromagnet) located at {..Xj..} is given by 

i^{xuX2,...,XN) = {-l)^^^X{&lTi^'\(t^{x,-X,)) (2) 

g; . 

I , with (3 — A. The model is interesting from several points of view, such as the connection with Gutzwiller projection 
'"O in strongly correlated systems, and from the intimate connection with the isotropic Heisenberg antiferromagnet, the 
g , Bethe chain. The spin- spin correlation function of the above wave-function at /3 = 4 has the same decay exponent as 
the Bethe chain, namely unity. Further, at /3 = 2, the wave function is the ground state of the free Fermi gas, with 
either long ranged hops or just nearest neighbor hops. The first case at /3 = 2 corresponds to dropping the zz part of 
Eq.(|l]), and the second to the anisotropic Heisenberg model: H = J^ti'^f'^f+i + ^f^f+i + AS'f Sf+i) at A = 0. Thus 
we find that /3 = 4 and 2 are in close correspondence with A = 1 and respectively. The Heisenberg model is well 
known to have a transition to a massive Neel ordered phase at A = 1, and so one might suspect that the wave 
function Eq.(||) could develop long ranged order as (3 increases from 4, perhaps even at /3 = 4+, a possibility we shall 
investigate in this paper. It is obvious that Neel order can also be viewed as crystalline order of hard core bosons, 
where the bosons correspond to the spin reversals of the Heisenberg system via the familiar lattice gas analogy. We 
will almost exclusively use this point of view below. Note that the density variables pi = (^ — S'f ) take values 0,1 so 
that we can map density correlators to spins readily, with N the number of hard core particles restricted to < L/2. 

For the case of /3 ^ 2,4 the wave functions do not represent either symmetric or antisymmetric functions, and are 
hard to interpret as physically allowed states for bosons/fermions, unless one imposes a rather non-analytic restriction 
of taking the modulus. For larger even integer values of /3 = 6,8.. the wave function Eq.(p|) is an eigenstate of the 



anisotropic version of the Hamiltonian Eq. (Ill) , but only in restricted sectors of numbers of 



^articles, for fillings up to 



^, since beyond this filling the states are no longer "good functions" in the sense of Ref. |^], i.e. they have Fourier 
components that "spill out" of the first Brillouin zone, requiring umklapp. 

The evaluation of correlations in the above wavefunction reduces to those of a 2-d Coulomb gas confined to a 1-d 
ring, but with the positions of the particles discretized to a lattice. This is a far reaching distinction from the case 
where the charged particles are in the continuum, a case that is familiar from the well known results of Dyson, Mehta 
and Gaudin for random matrices. In the latter case, the Coulomb gas does not crystallize in the sense of possessing 
LRO, although the density correlators are arbitrarily slowly decaying. In the lattice case one expects LRO, which is 
consistent with Neel order. 
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The discrete Coulomb gas has been subject to a few exact calculations earlier. Gaudin |^ computed the normaliza- 
tion constant of the wavefunction and the grand partition function exactly at three values of /? = 1, 2, 4. His isothermal 
calculation of the grand partition function at these values of the temperature gives the distribution of zeroes in the 
thermodynamic limit as lying on segments of the unit circle. Mehta and Mehta computed the density correlators 
exactly at these values of f3. The calculations at (3 = 2 are not unexpected, since the model at /3 — 2 reduces to a free 
Fermi lattice-gas, but the other cases are highly non trivial. Sutherland Q has presented results at zero temperature 
for the allowed ground state patterns, that turn out to be quite complex for arbitrary rational fillings. 

In this work, we present some asymptotically exact results on this problem, using a combination of renormalization 
group, and exact duality arguments on related models. We only consider simple rational fillings in this work with 
filling / = x = 1/2, 1/3, 1/4... We find that for each such filling f there is a transition temperature /3c = 2//^ at which 
the LRO sets in. In Section 2 of this paper, we show that the discrete Coulomb gas (DCG) is asymptotically equivalent 
to a sine Gordon model that has been studied extensively in connection with several interesting problems recently. 
This is achieved through a series of approximations that lead to a phonon representation at high temperatures, and 
a kink representation at low temperatures. Correlation functions are discussed in both representations. In Section 
3, the analogs of the phonon and kink representations are constructed for a model asymptotically equivalent to the 
discrete Coulomb gas, and an exact duality connecting the two pictures is obtained. The duality found by us is closely 
related to that found by Kjaer and Hilhorst who studied a discrete height problem- a roughening model, where 
the heights interact via a 1/r^ interaction. Our sine Gordon model reduces to this model on integrating out the 
Gaussian displacements. Section 4 presents numerical results that confirm the results of the previous sections and the 
difficulties involved in extracting true exponents are highlighted. Connections to related models and other general 
issues are discussed in Section 5. 



In order to understand the spin correlations for these wavefunctions, it is convenient to convert the problem into 
one in classical statistical mechanics. For a ring with N ~ f L spins, if ^{xi, xi . . . xn) is the amplitude for the down 
spins located at xi,X2 ■ ■ -Xn, then \ip{xi,X2, ■ . ■Xn)\^ is the corresponding probability. By rotational symmetry, the 
spin-spin correlation function {S{x) ■ S{y)) is equal to 3Sz(x)Sz{y)- This can be calculated from without any 
knowledge of the phase of ip. (This is not possible in general for higher order correlations, where one cannot always 
get rid of 5+ and S- operators by symmetry arguments. In this paper, we shall only consider two-point correlations.) 

If we express \tp{xi,X2 ■ ■ ■ xn)\'^ as exp[ln|V'P], we can view — Inj'f/'p as the energy of a classical system of N 
particles (distributed over L = N/f sites), and \tp\'^ as the statistical weight assigned in thermal equilibrium to a 
particle configuration. For the wavefunctions we consider here, — In |7/>p has the form 



The additive constant at the end of the right hand side is necessary in order to ensure that the wavefunction is 
normalized; however, in the statistical mechanics picture, it only gives rise to an overall proportionality factor in the 
partition function, and can be ignored. The l/c/)'^ inside the argument of the logarithm has been pulled out from the 
additive constant so that the L — > cx) limit exists. 

Eq.(||) describes a collection of particles with pairwise interactions. Every particle repels every other particle 
logarithmically. The argument of the logarithm is effectively (the square of) the straight line (or chord) distance 
between two points Xi and Xj on the ring. Thus the system consists of a charged gas with two-dimensional Coulomb 
interactions, but with the particles confined to a (one-dimensional) ring lattice. The prefactor /3 has a natural 
interpretation of the inverse temperature. 

In its ground state, the system has one particle on every 1//*'* site, at least for simple fractions / that we consider 
here ( say 1/2, 1/3 etc.). As the temperature is lowered, i.e. (3 is raised, there is the prospect of the system crystallizing 
into a long-range ordered state. As we shall see in this paper, this indeed happens at /3 = 2//^, and is the result of 
a combination of two factors. Firstly, althovigh one normally expects short range order in a one-dimensional system, 
the long ranged logarithmic interactions convert this to quasi long range order even at high temperatures. Secondly, 
the restriction that particles can only be placed on lattice sites crystallizes the system at low temperatures. 

There are two complementary approximations that one can make on the ring-lattice Coulomb gas, one appropriate 
for high temperatures and one for low temperatures. Both these approximations yield a one-dimensional long-ranged 
sine-Gordon model, but with different parameters, reflecting the well known duality of this model. This duality is 



II. PHONON AND KINK REPRESENTATIONS 
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usually derived for a continuum sine-Gordon model, Ijl^l but there are slight differences for the lattice version, as we 
shall now see. 

A. Phonon representation 

At high temperatures, there are large fluctuations in the separation between neighboring particles. It is reasonable 
to expect that, in this regime, the underlying lattice constraint might not be very important. Accordingly, in Eq.(^, 
we express Xi = + where Ui/ f is the deviation from an ideal crystalline state. The hard lattice constraint is 
replaced with a periodic potential, to obtain 



^2 fL 



The Ui's are now continuous variables, and the potential V favors locating the particles at lattice sites. As an 
example consider V{ui) = —\const\ cos(2Trxi) which reduces to —\const\cos{2'!TUi/f), leading to the sine Gordon 
theory considered below. More formally, the hard lattice constraint can be expressed in terms of a singular periodic 
potential V: V{u) = In exp(27rniu//)] . We shall instead consider the general class of potentials with periodicity 
/, and later exploit the universality under the renormalization group. Expanding the first term on the right hand side 
of Eq.(Q), to second order in u we have 

In Fourier space, expanding Uj = '^q{uq/\/N) exp[i j q], it can be showij^ that this is equivalent to 

q i 

The sum over q ranges from — vr to tt. Compared to the 7r|g| term, the q^ is irrelevant in the renormalization group 
sense, and can be ignored in calculating long wavelength properties. Higher order terms in the expansion of Eq. (^) in 
powers of u are similarly unimportant, as we see later. 

B. Coulomb Gas 

We now cast the problem of computing the partition function for Eq.(^) in the form of another Coulomb gas, but with 
a variable number of charge pairs, controlled by a chemical potential. Writing a general expansion exp ( — V{uj)^ = 
X^m ^"'h exp(27riujmj//), consistent with the periodicity u — > u + / we consider the partition function 

Z= /'Hexp[-^G(g)u,2_,]n(E^™.^'"™'''"'^0 

q j nij 

where the uj's are continuous variables, and G{q) is some function of q. (Corresponding to Eq.(^), one would have 
G{q) = (/3/2)(7r|g| - q^/2).) This is equivalent to 

^ = E E • ■ • / 1'^-] [ - E G{q)uqu-q] ( n -rr,y^' "^'"^ . (8) 

mi rn2 q j 

If G{q = 0) = 0, the only terms in this sum that are not zero are those for which J^j = 0- Integrating out the m's, 
we have 



^We use (f )2 EnJi^cosec2(^) exp(ifcn) = -lfc|7r + ifc^ + ^(1 - ^) with -tt < fc < tt. 
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(9) 



This is the partition function for a charged gas. (If we restrict Cm to be non-zero only for m = ±1 or 0, we have a 
dilute gas of unit positive and negative charges.) Using the fact that rrij = 0, the exponent in the exponential is 
equal to 
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1 - COS q{i-j) 



Giq) 



= 72 



dq 



1 - cosg(^ - j) 
Giq) 



(10) 



where we have taken the L — > oo limit in the last step. For the case G{q) = ^\q\ , i.e. the leading low energy part 
of Eq.(^, we can evaluate the integral easily for large separation of the charges, and find the dilute Coulomb gas 
partition function: 
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exp 



i<3 
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with j3ef f — 1^ and the 'chemical potential' /ie// — ~(7 + log7r)/2 = —.860973 with 7 Euler's constant (.577216). 

The Coulomb interaction binds unlike charges, and repels like charges. The quadratic terms of G{q) give rise to a 
sublogarithmic part in the interaction energy, but do not affect fJ-eff- The object fJ-eff is not really a chemical potential, 
since the relative weights assigned to different charges depends on the coefficients {cm}- However, for large Peff, where 
the partition function is dominated by rrii = 0, ±1, Heff may be viewed as the chemical potential associated with the 
(unit) positive and negative charges in the system, provided that c±i — cq. Alternatively, /ie// can be absorbed in a 
redefinition of the coefficients Cm, and is included when we start with general periodic potentials in Eq.(^. 



C. Kink representation 

In order to develop this representation, it is convenient to first rewrite Eq.(^ as 

const (12) 



f3H[{p}] J2 m^m, In 

{i<j} 



. 7i'(i - j) 
sm 



The variable nii is — f)/ f if the site i is occupied, and —1 if the site is unoccupied. Compared to Eq.(|^), this 
is equivalent to adding a background charge of — / at every lattice site, and then reducing the unit of charge to /; 
although this changes the total energy (the additive constant in Eq.(p^) is different from that in Eq.(||)), the energy 
difference between different configurations — and therefore their relative statistical weight — is unaltered. 

At low temperatures, the system is in an almost perfect crystalline state, with one particle after every 1/f sites. 
There are long segments that are shifted by I sites with respect to a reference crystalline state, with 1 = 0,1,2... (1//— 
1). (Shifting a segment by 1/f sites is equivalent to not shifting it at all.) There is an effective buildup of charge at 
the "kinks" at the segment boundaries. At low temperatures, the dominant configurations have neighboring segments 
with a relative shift ±1 corresponding to a unit positive charge (a kink) or a unit negative charge (an antikink) at 
the segment boundary. Note that there is no alternation rule for the kinks here, unlike the case for a system where 
the ground state has ferromagnetic order rather than antiferromagnetic. 

We now show that the energy of the system can be understood as an interaction of these charged kinks. As 
illustrated in Figure |^, two neighboring segments can be viewed as each consisting of a long string of quadrupoles, 
with a residual charge between them. Each quadrupole consists of / — 1 negative charges terminated by a charge 
+(/ ~ l)/2 at either end. The residual charge associated with the kink is ±1; in Figure |] , it is +1. Since the 
interaction energy between a quadrupole at i and a charge at j decays as l/(i — j)^, the total interaction energy 
between all the quadrupoles in a segment and a kink not at its terminus decays as 1// when the typical segment size 
is I. The interaction between the quadrupoles in non- adjacent segments can likewise be neglected. One is left with 
the interaction between the kinks. (The interaction between a kink and its adjacent segments, or two neighboring 
segments separated by a kink, can be interpreted as a self-energy or chemical potential for the kink; as in the previous 



section, we need not keep track of this.) The final picture that emerges is Eq.(ll), with rrii restricted to 0,±1 and 
Peff replaced by jSp. 
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FIG. 1. Typical lattice segment for / — 1/3. Every positively charged site is doubly charged. The vertical lines denote the 
boundaries between adjacent quadrupoles, and lie on the midpoints of the positively charged sites. The thick section in the 
middle of the line segment is a kink, consisting of one positive charge at either end and only one negative charge inside. The 
total charge of the kink is thus +1. 



D. Relationship between the representations 

The kink antikink gas, that we derived from Eq.(^ by a sequence of approximations, thus leads to a partition 
function that is essentially the same as obtained in the phonon representation, except that f3f'^ is replaced by 4//3/^. 
This is not surprising, since the continuum one-dimensional long-ranged sine-Gordon model has such a duality. pT| ] 
However, in the discrete version, when the strength of the potential V{ui) in Eq.lQ) tends to infinity, corresponding 
to a lattice model, the fugacity of the kinks in the kink representation does not tend to zero. This is because Ui can 
jump from one integer to another when the discrete variable i is increased by unity, whereas discontinuities in u{x) 
are not allowed in the continuum version. 

For (3 = 2//^, the system is at its self dual point. Of course, the duality that we have arrived at is only an 
approximate one. In particular, if one starts from Eq.(|^) at the self-dual (3 = 2//^, and proceeds along the phonon 
and kink routes, the resultant chemical potentials and various sublogarithmic interaction terms are different. In 
Section 3, we shall obtain an exact duality for a model using a restricted class of periodic potentials V{u), and a 
phonon spectrum that is linear in |g| only for small q. 



E. Renormalization Group 

The long- wavelength limit of Eq. (^) has been studied using the renormalization group [pT| in connection with several 
different problems. [ p^Jl3| We shall cite the results here without deriving them. 

From the \q\ form of the propagator for small q, one can see from power counting that the u{x) field is dimensionless in 
the absence of loop corrections. This is similar to the Kosterlitz-Thouless transition — or the sine-Gordon model — 
in two dimensional systems. Unlike the case there, however, the singular form of the propagator here prevents any 
renormalization of u (or equivalently of /?) to any loop order. 

If the potential V{u) is expanded in its harmonics, all higher harmonics die away rapidly compared to the lowest 
one, and can be neglected. Replacing V{u) with g cos(2ttu/ f), one obtains a one-dimensional long-range version of 
the sine-Gordon model. Calculating one-loop corrections, one finds dg/dl = g{l — 2/f3p) + 0[g^) (shifting u by 
//2 shows that there \s a, g ~* ~g symmetry) and hence the operator g is irrelevant for (3 < 2//^, and relevant for 
(3 > 2//^. Although this is a weak-coupling result, it has been argued to be true even for large g. At /? = 2//^, by 
mapping the problem to the scattering from a potential of one dimensional free fermions, it can be shown that the 
behavior is not universal, and depends on g. 

In the high temperature phase, Eq.(^ thus renormalizes to a harmonic phonon energy. The form of the density 
density correlations can be obtained from the following argument. The deviation of the density from the mean has a 
lowest Fourier component of the form Sp{x) = cos{2Trfx + 9{x)), where 9{x) varies more slowly than the oscillations in 
cos(27r/a;). Comparing with the definition of the (coarse-grained) displacement field u{x), we see that this is effectively 
cos [27r(/x -I- u(a:))] . The connected part of the density correlation function is then of the form 

Kc{x-y) ^ {5p{x)5p{y)) - l^cos 2tt [f x + u{x)) cos2ti [fy + u{y))y (13) 

Using the fact that u{x) is a Gaussian field, this simplifies to 
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K,{x -y)^ cos 2n{fx - fy) exp | - 27r2 ( [u{x) - u{y)] ' ) | ^ | ji^^|4//'^^ (14) 

where the last expression is an asymptotic resuh. Thus Kc{x—y) is a product of two terms: a rapidly oscillating factor, 
corresponding to the periodicity of the ideal crystalline state, and a factor that decays as a power of the separation 
between the points. In the low temperature phase, where g is relevant, u{x) is almost always at the minima of the 
potential V{u), and there is long range order in the system. This is easiest to see in the kink representation. 

For P > 2//^, proceeding through the kink representation yields Eq.(|6|), with /3 A/{Pf^), i.e. in the high 
temperature phase. As we shall now see, the resulting irrelevance of the operator g implies long range order in 
the kink representation. The relative phase of the local periodic structure at two points x and y, compared to a 
reference crystalline configuration, is 2Trf{nk{x,y) — n-i^{x,y)}, where nk{x,y) and n-j^{x,y) are the number of kinks 
and antikinks respectively between x and y. As in the phonon representation, the correlation function is then of the 
form 

Kc{x - y) ^ (^exp 2Trif{nk{x,y) - n-j:{x,y)} ^. (15) 
Since in the kink representation the potential V{u) in Eq.(0) with {3 — > 4//3/^ is interpreted as the generating function 



for the kinks and antikinks, the right hand side of Eq.(15) can be obtained by changing V{u) in the region between x 
and y from g cos{2ttu/ f) to gcos(27r/ + 27rw//) (so that p^'^'^^'^l S picks up a factor of e^"^^"^^). If the partition function 
thus modified is denoted by Z{g,g'), and the original partition function by Z{g,g), we have 

^exp 2mf{nk{x,y) - n-j:{x,y)] ^ = '^^{g g) ' ^^^^ 

If g is an irrelevant operator, this flows to a constant at long distances under renormalization (the value of the constant 
depends on the finite corrections that are removed along the course of the renormalization flow), and there is long 
range order. By expanding the right hand side of Eq.(^6|) in powers of g for small g, it can be seen that the correlation 
function decays to its long distance limit with a power law transient rather than an exponential. 

An alternative way to understand the correlations in the kink representation is to start with Eq.(|l^). At low 
temperatures, the system consists of bound pairs of kinks and antikinks. The probability that the number of kinks 
and antikinks between two points x and y separated by a large distance are not matched, is then dominated by cases 
when a kink (antikink) lies just inside the interval (x, y) and its partner lies just outside. This is clearly independent 
of the separation between x and y when the separation is large, so that K^ix — y) goes to a non-zero limit for large 
|a; — y\. The phase transition then corresponds to an unbinding point, where the mean separation between bound 

pairs diverges, i.e. /'"^ ^' dl (l/l^f^) diverges for large | X — y\. This occurs at f3 — 2/ p. Beyond this point, there is 
a proliferation of kinks and antikinks; for large separations mod(nfe — nj:, 1//) is equally likely to assume any of its 
allowed values, and there is no long range order. 

In more detail, the density density correlation function can be expressed in terms of the set of functions 

C^^Hj-k)^ {exp{27riiy{u,^uk)}), (17) 

with V = 1,2 ... . The preceding discussion only deals with = 1; other values of v give rise to corrections to the 
correlation function that decay more rapidly, and therefore do not affect the leading asymptotic behaviour. 



III. EXACT DUALITY IN VILLAIN-SINE-GORDON THEORY 



In this section we consider a particular type of sine Gordon Theory that corresponds to a Villain approximation 
of the cosine function, hence the Villain sine Gordon model (VsG). The advantage of this model is that one has 
an exact duality reminiscent of the Kramers Wannier duality in the 2-d Ising model. Towards this end we begin by 
considering a model for the energy in the sense of models described in Eq.(^, a sine Gordon model given by 

pH,G ^9^"? + ^9 [1 - cos(2^^/j//)] . (18) 

The Gaussian propagator hq is specified partly by giving its leading behavior as hq = \q\ + 0{q^) for small q, the sum 
is over the N wavevectors q, obeying — tt < q < it. The partition function is obtained by writing 
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ZsG= j IIj duj exp(-/3 Hsg) 
The Villain version of this model is defined by the partition function 

ZvsgIP, 9\= ^ I duj cxp - ^ XI ^1^-1 - 7^.9/3 ^ ^ ^j)^ 

c n_Li'' L r, 



(19) 



corresponding to a periodic function replacing the cosine in Eq.(18), as usual, with the correct quadratic coefRcient. 
The VsG model is defined by the above partition function, and the rules for computing the correlation functions given 
in Eq.(|r^), in terms of which the original density density correlation function of Eq.(||) can be expressed. (Formal 
expressions for C^i^J - for the VsG model are given in the Appendix, but the renormalization group arguments 
of the previous section are sufficient to obtain the qualitative behaviour.) The evaluation of the partition function is 
done by two different ways, leading to the same kink partition function, but with different parameters, and hence the 
duality follows. 

The first method is similar to the phonon representation described above, and is based on the Poisson summation 
formula: 



«j=0.±l,.. 



1 



E 



exp 



m,-=0,±l, 



2TTimjUj 
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. This is substituted in Eq.(|19D, leading to a shifted Gaussian in the variables Uq, which can be integrated out, yielding 
a product form for the partition function 



Zvsg[P,9] = (27rg/3) Zgauss{l3, {K}) Z^ortexiP, g], 



with 



^gauss{P,{hq}) = H J — 



2 

l3hg 



1/2 



(20) 



(21) 



With rUg defined as l/\/]V^ exp(«(7 j 



''j. 1 ^vortex 



is given by 



E 



m,=0,±l, 



exp 



2tt 



E 



(22) 



where the constraint ^ nij — 0, arises from the vanishing of hq at small q. 

The second method is as follows. In Eq.(|l9|) we fix a set of {^_,} and then shift the variables Uj — Uj + f 
Next, we integrate out the (still Gaussian) variables u'j using 



{u'qU'^q) 



1 



1 



j3TT hq + a 



(23) 



where a is defined as 



and find a factorization 



47r5' 



ZvsgW, g] = Zgauss [P, {hq + a ^}] Zr ough 



-'rough 



«j=0,±l,. 



(24) 



(25) 



The second part of the above is best seen as a roughening model with discrete height variables ~ 0,±1, ... interacting 
with a potential that is long ranged l/?'^, since the propagator is linear in \q\, for small \q\. In order to proceed, we 
define dual variables: rjj = ^j+i — ^j. These satisfy the condition J^Vj = ^ due to periodic boundary conditionss and 
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have a natural interpretation in terms of height differences in the roughening model. In terms of Fourier components 
iVqj S,q) = l/v^X) 6^P(*9 ' have rjq = {exp{iq) — 1) , and so the roughening model becomes 



Zrough \P , g\ 



?j=0,±l,.. 



(a hq + 1)(1 - cos(g)) 



(26) 



The propagator is now again of the form yiy for small j^l, and hence we can hope to get a more exact equality. This 
motivates us to choose the function hq — > h*q{a) which satisfies a quadratic equation: 



l)(l-cos(g)) 



4 



27ra 

7^ 



27r 



Thus, provided a < 4, we have a self dual propagator 



hq = 



2 


sin(g/2) 




1 - 


2a 


sin(g/2)| 



It is readily seen that h*q = \q\ + 0{q^) for small \q\. This choice gives 



Zrough[l3,g\^ ^ CXp 
«j=0,±l,.. 

4 

w 



27r 



VqV-q 



— z,, 



(27) 
(28) 

(29) 
(30) 



The second equation follows from comparing the first with the definition of the vortex partition function Eq.(^2|). For 
the self-dual propagator, ft,*, we see that Zrouoh and Zyortex are in fact independent of g, and are equivalent to the 
model studied earlier by Kjaer and Hilhorst. [g| However, the original ZvsG is still g dependent. The implications of 
this will be discussed further in the last section of this paper. Hence we have the final result: 



ZvsG[l3,g] = Zgauss[f3,{h* +a ^}] Z^ 



A comparison of the above with Eq.(|o|) provides the exact duality relation for this choice of h* 

1/2 



z„ 



H„ 



27rgf3hq 



(31) 



(32) 



There are several comments to be made at this point. Firstly the restriction a < ^ implies that the coupling 
constant g must be large enough in Eq.(19); too weak a periodic potential would have large fluctuations in Uq, that 
are not acceptable to this relation. In the limit of infinite g the relation is particularly simple, namely a = and 
hence h* — 2| sin(|)|. In this case the Villain approximation also would be exact for the sine-Gordon model Eq. 
The series of equivalences that have been established here can be summarized in the following diagram:. 



Zvsg[P,9] 
t 

Z rough 

7 4 
^vortex p j4 



T 

^rough 

t 



4 



(33) 



Zv 



sG 



0r 



■,9 



Here the symbols and < — > symbolize relations via Eqs.(p9|), ( |25| ) and ( pO| ) respectively. We see that the critical 
temperature of the model, if unique, is constrained to be /S^ — 2//^. There is of course no guarantee that there 
is no other critical point, if so, they must occur in pairs and satisfy the product condition PicP2c = (2//^)^. As 
mentioned earlier, it has been argued from renormalization group considerations that there is only a single critical 
point; numerical evidence is presented in the next section. 
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IV. NUMERICAL RESULTS 



The results in the previous sections have been based on perturbative calculations in the parameter g to obtain its 
relevance or irrelevance under renormalization. Although it has been argued that such perturbative considerations 
are in fact valid for all g [pl| for the continuum sine-Gordon theory, it is nevertheless useful to compare the results 
with numerical simulations, since the large g regime could be different for the continuum and discrete models. 

Correlation functions for the discrete Coulomb gas on a ring were computed numerically. Only the case of half- 
filling (/ = 5) was considered. The numerics were performed by starting in an ordered configuration, and evolving 
the system under Monte Carlo dynamics for some temperature. Various values of f3 were chosen, starting from /? = 4 
on the high temperature side, to /? = 9 on the low temperature side. This range covers the freezing transition at 
/3 = 2//2 = 8. 

Figure |^ shows a log-log plot of the connected part of the correlation function, Kc{r), as a function of r for (3 ~ A. 
For convenience of representation, a factor of (—1)'^ has been removed from Kc{r), so that it represents the deviation 
from perfect crystalline order. Lattice sizes of L = 30 through 960 were used. 
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FIG. 2. Log-log plot of the correlation function for /3 = 4. System sizes range from 30 to 960. Error bars are smaller than 
the symbols (except for the last three points). The dashed line corresponds to a power law with exponent —1, as per the exact 
resuh of Mehta. 

The correlation function is described very well by a power law decay with an exponent 1: Kc{r) ^ l/r. This is in 
agreement with the exact result. |^] 

For any /3 < 8, we expect a power law decay of Kc{r) with an exponent 4//3. However, Figure || shows a log-log 
plot of Kc{r) for P — 6; the apparent slope is significantly different from 4//3 — 2/3. 

One has to be careful in interpreting this result, since as /? — + 8, the irrelevant operator g renormalizes to zero more 
and more slowly, so that corrections to scaling could affect the apparent exponent over a fairly wide regime. To test 
whether the slope in Figure |^ can indeed be explained by leading irrelevant corrections, we try to fit the correlation 
fmiction to the scaling form Kc{r) ^ r^'^/^Kc{r/L]gr^^/^). This is the specific case for /? = 6 of the general scaling 
form 

Kc{r)-^k,(rmgr^^-^'^^), (34) 
based upon the RG flow of g, namely dg/dl = g{\ — 8//3) -|- 0{g^), with Kc{a, b) possessing a regular expansion for 
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FIG. 3. Log-log plot of the correlation function for (3 — 6. System sizes range from 30 to 960. Error bars are smaller than 
the symbols. The dashed line corresponds to a power law with the best fit exponent —0.553, which is significantly different 
from the theoretical result. This can be understood in terms of corrections to the scaling form, as shown in the next figure. 

a,b ^ 0. If Kc{r)r^^'^ is obtained for r = XL for fixed A and varying L, the result should then be a function of A and 
gj,-i/3_ small g, one would expect this to be a linear function of r^^/"^. Figure |^ shows such a plot oi Kc{r)r^^^ as 
a function of r~^/^ for various values of A. A set of straight lines is obtained, consistent with the scaling prediction. 
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FIG. 4. Correlation function for /3 = 6, multiplied by its asymptotic power law decay. The x axis plots the r-dependence 
of the leading irrelevant operator, g. The different symbols correspond to r = XL with different values of A. Apart from the 
smallest values of r, the data fits reasonably well to a set of straight lines. Note that r decreases along the x axis. 

In view of the strong dependence on /3 of the dimension of the operator g, it is not very useful to try this for larger 
values of /?, since the range one obtains for L^~^^^ is quite limited. Conversely, there is no sign of any corrections to 
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scaling for /? = 4, because of the rapid decay of g under renormalization. 

At /? = 8, the operator g is marginal. This is not just a perturbative result; based on the continuum model, one 
expects a line of fixed points for /3 = 8, with continuously varying asymptotic behavior as the initial g is varied. In 
order to check this scenario, we performed Monte Carlo simulations on a slightly modified version of the DCG at half 
filling. The lattice size was doubled, corresponding to / = |, but the particles were biased to be on the even sites of 
the lattice by adding an extra potential to the odd sites. It is clear that if the bias is infinite, the system is equivalent 
to the DCG at half filling, while with zero bias, one has the DCG at quarter filling. In general, the lattice structure 
can be represented as a strong (strictly speaking, singular) potential W{u) that is periodic under m ^ u + |, and an 
additional weak potential V{u) that is periodic under m — > u + i. The strength of V[u) depends on the bias favoring 
the even sites. Since W{u) is irrelevant for 13 < 32, only the potential V{u) affects the asymptotics. As one adjusts 
the bias, which corresponds to changing g, one should see a continuous evolution in the asymptotic behavior of Kc{r). 
Figure || shows that this is indeed the case for L = 960. 
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FIG. 5. Log-log plot of the correlation function for f3 — 8, and system size L = 960 showing the effect of varying the initial g, 
a marginal operator. Error bars are smaller than the symbols. The different curves correspond to different values of the bias, 
which is the relative 'Boltzmann' weight of the odd sites compared to the even sites. The correlation function is computed 
by first coarse graining the density, so that there are effectively 480 sites (corresponding to 240 particles at half filling), and 
removing the factor of (—1)*^. As the bias changes, the correlation function evolves smoothly, with no universality seen even for 
large r. 

Beyond the transition. Figure H shows the correlation function for (3 = 9, indicating that long-range order has set 
in. The slow decay of the operator g in the dual representation (since 64/(3 is not much less than 8) leads to the long 
transients in the correlation function. 
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FIG. 6. Log-log plot of the correlation function for (3 — 9. System sizes range from 30 to 960. Error bars are smaller than 
the symbols. The dashed line is of the form A + Br^^^^ , which includes the leading scaling correction to the long range order. 
The parameters A and B are adjusted for approximately the best fit to the eye, and are both equal to 0.1. 

By comparison, the correlation function at P — 10, shown in Figure approaches the asymptotic r —> oo limit 
much more rapidly. 
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FIG. 7. Log-log plot of the correlation function for /3 = 10. System sizes range from 30 to 960. The dashed line is of the 
form A + Br"^^^, which includes the leading scaling correction to the long range order. The parameters A and B are adjusted 
for approximately the best fit to the eye, and are chosen to be A = 0.189 and B = 0.032. With these parameters, C(cx))/C(0) 
is estimated to be approximately 0.75. 

We did not increase (3 beyond this, because we do not expect to see any quahtative change in the correlation 
function, and because as /3 is increased, it takes progressively longer for the system to equilibrate. 
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Thus we see that the numerical results at half fiUing for the density density correlation function agree on both sides 
of the freezing transition with the analytical results obtained in the previous sections. 



V. DISCUSSION 



In this paper, we have obtained the phase diagram for the two-dimensional Coulomb gas on a one-dimensional lattice, 
which emerges as a generalization of the wavefunction of a spin ^ Heisenberg antiferromagnet in one dimension. We 
have seen that the gas freezes into a state with long range order below a freezing temperature Tc = P /2 that decreases 
as the density is reduced. Above this temperature, the gas has quasi long range order, with continuously varying 
exponents. The dependence of the LRO on /3 is thus quite dissimilar to the dependence of LRO upon A in the XX Z 
model. In the latter, the isotropic point is at the brink of crystallization, whereas in this problem, the wavefunction 
at /3 = 4 for the Haldane Shastry model ||l|,|| is well inside the power law phase, since crystallization only sets in for 
(3 > 2//^. On the other hand, the wavefunctions require "umklapp" for /3 > 2// beyond which the kinetic energy 
ceases to act in a simple way on these functions Q]. 

The strong connections that exist to earlier work |]9|- [l^ , |l6|| have been alluded to; we return to them here in more 
detail. In Section 3, we considered a one parameter family of Villain sine-Gordon models at each temperature, to 
obtain an exact duality transformation; the family was parametrized by g. Rather surprisingly, all the models in 
the one parameter family map to the same roughening model, independent of g. As shown in the Appendix, formal 
expressions for the correlation functions of the VsG model can be obtained in terms of the roughening model, with 
the parameter g affecting the expressions only at short distances. This implies that, for any inverse temperature /?, all 
the VsG models (independent of g) flow to the same fixed point under renormalization. Although this is reasonable 
at any other temperature, it is somewhat unexpected for /? = 2//^, where one has a line of fixed points. Although g is 
a marginal operator here, the self-dual propagator h*, and therefore the strength of the irrelevant operators, depend 
on g. One must conclude that, within the one parameter family that we construct, this change in the irrelevant 
operators is just sufficient to drive g to the same fixed point under renormalization, regardless of its bare value. Of 
course, one could obtain a duality mapping between VsG models with different, 'conjugate', choices of hq, chosen to 
satisfy 2(1 — cosg)(Q! + l/hq) — h'q/{l + a'h'q). However, although this would allow one to access other fixed points at 
P = 2//^, one would not have (strict) self-duality. 

As mentioned earlier, the roughening model constructed in Section 3 has been studied earlier by Kjaer and Hilhorst. 
1^ They obtained the phase diagram we have here, and at (in our notation) /3 = 2//^, exploit the self- duality to 
calculate the {mqun^q) correlator. As shown in the Appendix, the correlation function for the VsG model (and 
ultimately the Coulomb gas) is related to (expi[m(x) — m(?;)]). It is tempting to conjecture that m{x) can be treated 
as a Gaussian variable for the long distance form of this correlation function, especially because using the result 
of Kjaer and Hilhorst then yields a power law decay for t he V sG co rrel ation at the self-dual point with the same 
exponent, 1/4, with two different approaches given in Eqs.(A2) and (A3). However, numerical simulations we have 
conducted for the VsG model (/ = i,/3 = 8) do not bear this out: the numerical exponent is approximately 1/8. 

The continuum version of the long range 1-d sine-Gordon has been studied extensively in connection with dissipative 
quantum mechanics, ||l^ Luttinger liquids , ITT] and the quantum Hall effect. |l^ However, the duality transformation 
is slightly different from the discrete case. |lC|] For large g, u{x) must be close to an integer everywhere; a kink consists 
of a rapid change of u{x) from one integer to another. If (without a lattice cutoff) the theory is regulated with an 
mq^ term in the propagator for u, it is clear that u{x) cannot change discontinuously. The competing effects of g and 
TO then yield an effective kink size of ~ y/rn/g, and a kink fugacity of ^ exp(— /3^mg). This is equivalent to a sine 
Gordon theory expanded in powers of the cosine interaction, provided one chooses go ~ exp(— /Sy^mg). Thus we see 
that the large g regime maps to small g under duality. This is in contrast to the discrete model, where Ui can jump 
from one integer to another as i increases by 1, without any extra (i.e. not accounted for by hq) energy associated 
with the kink. Indeed, in the Villain version studied in Section 3, a large g maps to the same (large) g under duality. 

The renormalization group flows for the continuum model were obtained | [Tl| ] in the small coupling constant regime. 
Exploiting the duality transformation, it was possible to obtain the RG flows for very large coupling constants as well. 
It was argued that the RG flows could be connected smoothly between these two extremities; this was strengthened 
by showing that at the self-dual 'temperature', which corresponds to the scattering of non-interacting fermions from 
a barrier in the Luttinger liquid version, one should indeed have a line of flxed points. 

While for a small coupling constant the discrete and continuum versions should not differ in any physical way, 
this is not necessary when the coupling constant is large. Thus the possibility of a non-trivial strong coupling phase 
cannot be ruled out for the discrete model based on continuum arguments. Our numerics indicate that the Coulomb 
gas shows the same behaviour as the continuum model, suggesting that it is unlikely that there is such a strong 
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coupling phase. Notice that, although the RG flows are oppositely oriented for [3 > 2//^ and P < 2//^, under the 
duality transformation of Section 3 a larger coupling constant maps to a larger coupling constant, emphasizing the 
importance of the simultaneous change in the irrelevant operators. 

A one dimensional Coulomb gas with logarithmic interactions was considered earlier by Anderson, Yuval and 
Hamman [ pT[ in their study of the Kondo problem. This is equivalent to the kinks in a ferromagnetic Ising spin 
chain with long range coupling ~ 1/r^. However, the ferromagnetic nature of the underlying order forces charges to 
alternate. As shown by a real space RG calculation, jlj] integrating out tightly bound charge pairs renormalizes the 
strength of the logarithmic interaction between the remaining charges. Thus (3 flows under the RG, and one obtains 
the two-dimensional Kosterlitz-Thouless phase diagram. For finite coupling constant (or charge fugacity), this 
prevents one from obtaining the phase transition point exactly. 
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APPENDIX A: CORRELATION FUNCTIONS 



In this appendix, we obtain expressions for the correlation function C'^''(j) within the Villain sine Gordon theory 
using the two methods described in Section 3: the low and the high temperature limits. We begin by writing the 
correlation function explicitly as follows: 



C{k- l\(3,g) = (exp27ri(ufe - ui)) 



1 



ZvsGW,g) 



Pn/2^hqUqU-q - ^T2'^9^{U] - f + 27ri(ufc - Ul). 



(Al) 



f 



We can proceed to evaluate this by the two methods discussed above. 

The First method: We use the Poisson formula to trade the integer valued variables £,j in favor of the m'jS and 
find 



Cik-l;P,g) 



1 



ZvsG{P,g) 



E 



n," duj exp 



Ptt/2 hqUqU-q + 2-^1 / f m'j\ 



where m'^ ~ rrij if j ^ k,l and m'j. = rrik + f, m[ = mi — f. Now it is straightforward to integrate out the Gaussian 

^^(exp(igfc) — exp{iql)), we find 



3 

variables Uj and in terms of the variables fhq defined earlier and Srhq 



C{k - l;P,g) = exp 



exp - 



An 



E jr^^ ^ cos(gfc - ql)] 



(A2) 



47r 



y- 



W — 



niq Sm^q 



vortex[0,g] 



The prefactor decays as a power law at all (3, and the average in the second term is over the vortex partition function. 

The second method: We fix the variables £,j in Eq.( |Al| ) and shift the variables Uj — u'^ + f ^j, and integrate over 
Uj. We next use the difference variables rjj = — as in derivation of the roughening model Eq.(p^) , and find 
after some manipulations 



C(k - l;l3,g) = exp 



47r 



2m 



exp 



— {l — cos{qk — ql)] 



(A3) 



(1 + ahq){eyip{~iq) - 1) 



vortex[4/(/3/4),g] 
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The remarkable identity of Eq. ( [A^ ) and Eq. ( [A^ ) is a consequence of the two representations of the partition function. 
These are in general very hard to evaluate, the only simple situation is the case of very low temperatures, where we 
can assume a dilute gas of vortex anti vortex pairs. 
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